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ABSTRACT 

We present the stellar and substellar mass function of the open cluster 
IC 2391, plus its radial dependence, and use this to put constraints on the for- 
mation mechanism of brown dwarfs. Our multiband optical and infrared pho- 
tometric survey with spectroscopic follow-up covers 11 square degrees, making 
it the largest survey of this cluster to date. We observe a radial variation in 
the mass function over the range 0.072 to 0.3 M©, but no significant variation 
in the mass function below the substellar boundary at the three cluster radius 
intervals analyzed. This lack of radial variation for low masses is what we would 
expect with the ejection scenario for brown dwarf formation, although consider- 
ing that IC 2391 has an age about three times older than its crossing time, we 
expect that brown dwarfs with a velocity greater than the escape velocity have 
already escaped the cluster. Alternatively, the variation in the mass function of 
the stellar objects could be an indication that they have undergone mass seg- 
regation via dynamical evolution. We also observe a significant variation across 
the cluster in the colour of the (background) field star locus in colour-magnitude 
diagrams and conclude that this is due to variable background extinction in the 
Galactic plane. From our preliminary spectroscopic follow-up to confirm brown 
dwarf status and cluster membership, we find that all candidates are M dwarfs 
(in either the field or the cluster), demonstrating the efficiency of our photo- 
metric selection method in avoiding contaminants (e.g. red giants). About half 
of our photometric candidates for which we have spectra are spectroscopically- 
confirmed as cluster members; two are new spectroscopically-confirmed brown 
dwarf members of IC 2391. 

Subject headings: stars: low-mass, brown dwarfs, mass function — open clusters: 
individual (IC 2391) 
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INTRODUCTION 



The origin and evolution of brown dwarfs (BD) remains a fundamental open question. 
BDs have masses bridging the lowest mass hydrogen-burning stars and giant planets, so any 
picture of star and planet formation is incomplete if it cannot account for BDs. Several for- 
mation mechanisms have been proposed, including star-like formation from the compression 
and fragmentation of a dense molecular cloud, planet-like formation in a circumstellar disk, 
and the dynamical interruption of a star-like accretion process. 

There are observational signatures which may be used to distinguish between these 
scenarios, such as the distribution of binaries, the presence a nd properties of circ umstellar 
disks, the (initial) mass function (MF) and kinematics (see iLuhman et al.l 120071 for a re- 
view of observational signatures on the formation of BDs). Work over the past ten years 
has seen considerabl e success in measuring the MF into the BD reginie in several clusters , 



including a Orionis (iGonsalez-Garcfa et al 



the Orion Nebula Cluster (ONC) (IHillenbrand et a.1. 



(iBarrado v Navascues et al.ll2004l : iPlatais et al 



20061: 



Caballero et al. 



200G 



2007 



Lodieu et al.l 120091) . 



Slesnick et al 



20ol - IC 2391 



2007h and the Pleiades dMoraux et aP 12003 



Lodieu et al.l 120071 ) . These are only examples among many other analysis of stellar and 
substellar populations. 

The comparison of the MF in clusters with different p roperties (e.g. the different density 
clusters Taurus and ONC; clusters with different ages, e.g. IChabrierll2003l ) has led some work- 
ers t o draw conclusions abou t the relative efficiency of p ossible BD formation mechanisms 
(e.g. iKroupa &: Bouvierl 120031). While some observations (ILuhman et al.l 120071 ) and theoret- 



ical works (jPadoan fc Nordlundl l2004l : iHennebelle fc Chabrierl |2008| ) conclude in a common 
formation mechanism for BDs a nd stars, some studies has sugg; ested that BDs could form 
by massive disc fra gmentation fStamatellos fc WhitworthI 120081 ) . photoevaporati on of the 
accre t ion envelope (IHester et al.lll996l ). or inte rruption of the a ccreti on process (IReipurth 
2000l : iReipurth fc Clarkd l200ll ). For instance, iBate fc Bonnelll (120051 ) have performed hy- 
drodynamical simulations of star formation from fragmentation of molecular clouds. They 
concluded that objects which end up as BDs stop accreting before they reach the hydrogen 
burning limit because they are ejected from the dense gas soon after their formation by 
dynamical interaction in unstable multiple systems. 

This ejection scenario in some cases predicts a higher velocity dispersion and spatial 
spread of BDs in comp arison to stellar objects, w hich in turn may be observed as a variation 
in the MF with radius (IKroupa fc Bouvierll2003l ). On the other hand, other work have shown 
that i f BDs are fo rr aed by ejection , the v elocity distribution could be the same for BDs and 
stars (jBatel 120091 ). iMuench et al.l (120031 ) observed a radial variation in the MF of IC 348 
measured over 0.5 to 0.08 Mq, but no variation was detected in the BD regime. In a study 
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of the spatial di stribution of subste l lar ob jects in IC 348 as well as Trapezium in the Orion 
Nebula Cluster, iKumar &: Schmeja observed the stellar objects to be more clustered 

than the substellar ones, which they took as evidence in favour of the ejection scenario. By 



looking at the spatial distribution of the Taurus stellar and substellar population. iGuieu et al. 



( 120061 ) observed a gradient in the BD abundance relative to stars, which they conclude as in 
favour of the ejection scenario. 

In this paper, we present the results of a program to study, in detail, the MF of one of 
the nearest and ri chest open clusters, IC 2391. Th is has an age of 50Myr measured from 
l ithium depletion (Barrado y Navascues et al.ll2004l ) or 40M yr from main-sequen ce fitting 



( iPlatais et al.l 120071 ). The Hipparcos distance is 146.0_4;5pc ( iRobichon et al 



metalicity and extinction are [Fe/H]= —0.03 ± 0.07 and E{B — V) = 0.01, (iRandich et al. 



1999) and the 



2001). Because of its proximity and youth , this cluster has been subject to several studies 



Barrado V Navascues et al 



2004 



2001 



1999 



Koen fc Ishihara 



2006: ISiegler et al.ll2007 



Platais et al.l 120071 ). iBarrado v Navascues et al.l (l2004f ) and IPoddl tOQ4 \ measured the MF 
down to the substellar limit, but the MF for confirmed member with a completeness limit 
at lower masses has not yet been determined. 



Since IC 2391 is not as yo ung as IC 348 (age ~2Myr from iMuench et al.l l2003l ) and 
Trapezium (age~0.8Myr from iMuench et al.ll2002l ). one may expect it has already lost a 
significant fraction of its substellar population to evaporation by dynarn ical evolution. Using 
the tidal radius and mass of IC 2391 estimated by lPiskunov et al.l ( 120071 ) (7.4 pc and 175 Mq), 
we compute the escape velocity as Ve =0.4km/s and the crossing as time tcross =17Myr. (We 
stress that the crossing time is just an ord er-of-magnitude quan tity, 2R/v, where we have 
adopted for R the t idal radius of 7.4pc from lPiskunov et al.ll2007l and the velocity dispersion 
of ~0.85 km/s from lPlatais et al. I l2007h. Assuming a minimum value for the number of cluster 
members as the objects reported by lDoddl (120041 ) and lBarrado y Navascues et al.l (120041 ) (125 
and 33 respectively, together a total lower limit of 158 objects), we estimate that the lower 
limit of the relaxation time for this cluster is tj-eiax ~105Myr. F urthermore, in a nu merical 
simulation of open clusters and the population of BDs members, lAdams et al.l (|2002l ) shows 
that there would still be more than 80% of the original BD population in the cluster even 
after about 10 crossing times, assuming that BDs and stars have a similar velocity dispersion. 
Therefore, IC 2391 is still young enough for a radial study of its very low mass star (VLMS) 
and BD populations, consi dering the fact that ma ss segregation occurs on a timescale of 



order one relaxation time (IBonnell fc Davied 1 19981 : although recent work by lAUison et al. 



20091 suggest that mass segregation can occur on a smaller time than a relaxation time, at 
least for more massive stars). 



The paper is structured as follows. We will first present the data set, reduction procedure 
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and calibration in ^ We then discuss our candidate selection procedure in ^ present the 
survey results in §1] and then discuss the radial variation of the MF in §51 The preliminary 
spectroscopic follow-up is presented in §21 followed by our conclusions in ^ 



2. OBSERVATIONS, DATA REDUCTIONS AND CALIBRATIONS 

2.1. Observations 

The survey consists of 35 34x33 arcmin fields extending to 3 degrees from the center 
of the cluster and centered on RA=08:40:36 DEC=-53:02:00 (Figure [I]). The central 4 
fields will be referred to as the deep fields while the other 31 other fields will be referred 
as the radial fields and outward fields. (We make a distinction since they were observed 
with different exposure times and different filters, as will be specified below). The fields 
were chosen to extend preferentially along lines of constant Galactic latitude in an attempt 
to reduce systematic errors in any established cluster MF gradient which could arise from 
contamination by a Galactic disk population gradient. The total coverage of our survey is 



10.9 sq. deg. This compares to 2.5 deg^ in the survey of iBarrado y Navascues et al.l (120011 ). 



The optical observations were carried out in four runs with the Wide Field Image (WFI) 



on the 2.2m telescope at La Silla (iBaade et al.lll999l ) in : 24 January - 9 February 1999, 20 



- 24 January 2000, 10 - 23 March 2007 and 15 - 18 May 2007. The WFI is a mosaic camera 
comprising 4x2 CCDs each with 2kx4k pixels delivering a total field of view of 34x33 arcmin 
at 0.238 arcsec per pixel. The deep fields in our survey were observed in four medium bands 
filters, namely 770/19, 815/20, 856/14 and 914/27 (where the filter name notation is central 
wavelength on the full width at half maximum, FWHM, in nm) and one broad band filter, 
Rc. The radial fields were observed in Rc, 815/20 and 914/27 while the five outward fields 
were not observed in R^. 

These filters were chosen to sample the spectra of late M and early L dwarfs to improve 
selection over, say, Rc — Ic, and to minimize the Earth-sky background plus any nebular 
emission (as the filters are in regions of low emission). The pass band function for all 
filters are shown in Figure [2J For all radial fields, we have used an exposure time of 15, 10 
and 25 minutes for the 815/20, 914/27 and -Rc filters respectively. For each deep field we 
obtained an integration time of 65 minutes for each of 770/19 and -Rc and from 50 to 155 
minutes for 815/20, from 25 to 261 minutes for 856/14, and 50 to 80 minutes for 914/27. We 
additionally obtained short exposures for all fields to extend the dynamic range to brighter 
objects. The photometry from the short exposures were combined with the photometry of 
our long exposures for our analysis. In order to improve the determination of their low mass 
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status (via a better determination of spectral type and luminosity), we also observed all 
radial fields, including the outward fields, in the J-band using the Camera PAnoramique 
Proche-InfraRouge (CPAPIR) on the 1.5m telescope at Cerro Tololo, Chile (runs on 28 
February - 3 March 2007 and 10 March 2007). However, we did not get J-band photometry 
for the deep fields. CPAPIR consists of one Hawaii II detector of 2kx2k pixels for a field 
of view of 35x35 arcmin with a pixel scale of 1.03 arcsec per pixel. All fields were observed 
with a total exposure time of 30 minutes. The J filter of CPAPIR is centered at 1 250 nm 
with a FWHM of 160 nm. A detailed list of the fields observed with pointing, filter used, 
exposure time and lOa detection limit, is given in Table [H 

Our lOo" detection limit is J=17.7 and 914/27=20.5 for the radial and deep fields re- 
spectively (which corresponds to ~O.O3M0 for both cases). However, we can't expect to 
detect all objects down to these magnitudes. The completeness spanning from the brightest 
objects without saturation down to the lOo" detection limit is estimated by taking the ratio 
of the number of objects detected to the predicted number of detections and assuming a 
uniform distribution of stars along the line of sight in the fields of our survey. The predicted 
number of detections is derived very simply by extrapolating to the detection limit a linear 
fit to the histogram of the number of detections as a function of magnitude (Figure [3]). The 
completeness of the radial part of the survey is 91.8% while for the deep part it is 82.7%. 

The spectroscopic observations were carried out with HYDRA, a multi-object, fiber- 
fed spectrograph on the 4m telescope at Cerro Tololo on the nights of 6 and 7 January 
2007. Only two fields could be observed: the deep fields 15 and 20 (3 and 2 exposures of 
45 minutes respectively). We used the red fiber cable with the KPGLF grating (632 lines 
mm~^) and a blaze angle of 14.7° (no blocking filter was used) and . This gives us a coverage 
of 6429-8760 A centered at 7593 A and a spectral resolution of 4.0 A. 

2.2. Reduction and Astrometry 

The standard CCD reduction steps (overscan subtraction, trimming and fiat-fielding 
for the WFI data and dark subtraction, flat-fielding for CPAPIR data) were performed on 
a nightly basis using the ccdred package under IRAF. For WFI data we used the dome 
flat both for pixel-to-pixel variation correction and to correct the global illumination, while 
for CPAPIR data we used superflat (obtained by combining science image frames for each 
nights). For WFI data, we reduced each of the 8 CCDs in the mosaic independently and 
in the final step scaled them to a common flux response level. We made an initial sky 
subtraction via a low-order fit to the optical data, and for the infrared data by subtracting 
a median combination of all (unregistered) images of the science frames. Images were fringe 
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subtracted when fringes were vi sible, which was th e case for all medium bands filters used, 
in a similar way as described by Bailer- Jones et al. ( 2001 1^ Finally, the individual images of 



a given field were registered and median combined. We calculated magnitudes via aperture 



photo metry together with an aperture correction following the technique used in iHowell 



f 



( 1l989l ). An astrometric solution was achieved using the IRAF package imcoords and the 
tasks ccxymatch, ccmap and cctran. For each field, this solution was computed for the 
815/20 band image using the Two Micron All Sky Survey (2MASS) catalogue as a reference. 
The RMS accuracy of our astrometric solution is within 0.15-0.20 arcsec for WFI data and 
within 0.3-0.4 arcsec for CPAPIR data. For WFI data, the astrometry was also performed 
on a CCD-by-CCD basis. 

We only retained images for this study which were taken under photometric conditions, 
as determined by our monitoring of conditions during the observations and, moreover, by 
our data reduction procedure. 



2.3. Photometric Calibration 

To correct for Earth-atmospheric absorption on the photometry, we solved by least 
squares fitting the equation, 

rriA = rriA-raw + Za + CAiruA - ruB) + ^aXa, (1) 

for observations of the standard stars at a range of airmasses, where the spectrophotometric 
standard stars observed were Hiltner 600, HR 3454 and LTT 3864. (In order to obtain 
the o bserved magnitude from equations [H the fluxes of each standard star, /a, were taken 



from lHamuy et al.lll992l . 11994 ) The parameters and are the apparent magnitudes 
of our spectrophotometric standard in two particular bands {A and 5), where mA^raw is 
the instrumental magnitude of our spectrophotometric standard stars, Za is the zero point 
offset, Ca is the colour correction and the extinction coefficient for band A and X is the 
airmass at which rriA^raw was obtained. 

We calibrated the infrared data using the J band values of 2MASS objects which were 
observed in the science fields. By determining a constant offset between the magnitude of 
2MASS and our instrumental magnitude, we obtained the zero point offset. Since this zero 



"'^A fringe correction frame was created, which is a median combination of all science in a same filter with 
same exposure time. This fringe correction frame was scaled by a factor - determined manually for each 
science frames - and subtracted from the science image. 
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point offset was obtained with objects in the same field of view in each science frame, we 
did not perform a colour or airmass correction when reducing our NIR photometry. 



2.4. Mass and Effective Temperature Based on Photometry 

We used the spectral energy distribution to derive the mass and effective temperature, 
Tgff, assumi ng that all our photo metric candidates belong to I C 2391. We used e volutio nary 
tracks from lBaraffe et al.l (119981 ) and atmosphere models from iHauschildt et al.l (Il999al ) (as- 
suming a dust-free atmosphere; the NextGen model) to compute an isochrone for IC 2391 
using an age of 50 Myr, distance of 146 pc, a solar metalicity and neglecting the reddening 
{E{B — V) = 0). These models and assumptions provide us with a prediction of fx, the 
spectral energy distribution received at the Earth (in erg cm~^ s~^ A~^) from the source. 
We need to convert these spectra to magnitudes in the filters we used. Denoting as Sa{^) 
the (known) total transmission function of filter A (including the CCD quantum efficiency 
and assuming telescope and instrumental throughput is flat), then the flux measured in the 
filter is 



fA 



(2) 



The corresponding magnitude rriA in the Johnson photometric system is given by 



ruA = -2.5 log /a + ca, (3) 

where ca is a constant (zero point) that remains to be determined in order to put the model- 
predicted magnitude onto the Johnson system. We determine this constant for each of the 
bands Rc, 770/19, 815/20, 856/14, 914/27 and J in the standard way, namely by requiring 
that the spectr um of Vega produc e a magnitude, niA, of zero in all bands. Using the Vega 



spectrum from IColina et al.l (119921 ) we derive values of crc = —21.6607 C770/19 = —22.2517, 
C815/20 = —22.4391, C856/14 = —22.6341, C914/27 = —22.8353 and cj = —23.6865 mag. Apply- 
ing the two equations above to a whole set of model spectra produces a theoretical isochrone 
in colour-magnitude space. Note that this procedure only provides us with the "true" 
magnitudes of the model spectra, not their instrumental ones. The photometric calibration 
procedure applied to the data converts the measured, instrumental magnitudes to the "true" 
magnitude plane where we can compare them. 



Assuming that all our photometric candidates belong to IC 2391, we derive masses and 
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effective temperatures in the following way. We first normalize the measured spectral energy 
distribution (multiband photometry) of each object to the energy distribution of the model 
using the 815/20 filter. We then estimate the mass and effective temperature (which are 
not independent of course) via a least squares fit between the measured spectral energy 
distribution and the model spectral energy distribution from the isochrone. 

There are several sources of error in the mass and T^s estimates. These come from 
the photometry, the photometric calibration, the least squares fitting (imperfect model) and 
the uncertainties on the age of IC 2391 (we use 5Myr). This last and most significant 
error gives 0.075±0.006 M0 and 2914±43K for an object at the hydrogen burning limit and 
1. 000±0. 027 Mq and 5270±70K for a solar type object. 



2.5. Spectroscopic Data Reduction and Calibration 

The standard CCD reductions (overscan subtraction and trimming) were performed on 
each image using the ccdred package under IRAF. We then used the IRAF package dohydra 
to perform fiat-fielding (using dome fiats), throughput correction (with the skyfiats) and 
scattered light corrections. The spectra were wavelength calibrated using the PENRAY 
comparison lamp with 2 sec exposure time. Sky subtraction was performed in a similar 
manner as fringe subtraction in photometry: a standard sky spectrum (shown at Figure H]) 
was obtained from the median of our sky spectra (more than 20 fibers were assigned for 
sky subtraction in each Hydra pointing) and scaled to optimize the sky subtraction for each 
science spectrum individually. However, this sometimes resulted in Ha apparently being in 
absorption for some objects. We attribute this to Ha emission from the background itself. 
This is spatially variable and so subtracting the sky spectrum (which includes Ha) sometimes 
results in a over-subtraction of this feature. We discuss this contamination problem and the 
danger of determining membership status based on Ha in §6.11 

Finally, flux calibration was performed with the spectrophotometric standard Hilt- 
ner 600, which was observed three times a night, at three different airmass. 



2.6. Spectral type, effective temperature and mass determination 



For the objects for which spectra is av a ilable , we estimated in addition the spectral 
type using the PCS index from lMartfn et al.l (119991 ). The distinction betw een M-dwar f and 
background M-giants and M-supergiants was achieved using a CaH index (|Joneslll973l ). We 
visually inspected all spectra in order to conflrm the spectral type and luminosity class 
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estimati on. We e s timat ed a spectroscopic Teg from the spectral typ e using the te r apera ture 
scale of iLuhmanI (119991 ) for objects between MIV to M9V and of iMartfn et al.l (119991 ) for 
objects from LOV and later. We then use our isochrone of IC 2391 to obtain the mass based 
on Teff. 



3. CANDIDATE SELECTION PROCEDURE 

The selection procedure discuss here concerns only our photometric data while the 
discussion of the selection of our spectroscopic candidates is done in §61 The candidate 
selection procedure is as follows (and explained in more detail in the remainder of this 
section). Candidates were first selected based on colour- magnitude diagrams (CMDs). A 
second selection was performed using colour-colour diagrams. Third, astrometry was used 
to remove objects with high proper motion. Finally, non-candidates were rejected based on 
a discrepancy between the observed magnitude in 815/20 and the magnitude in this band 
computed with the isochrone of IC 2391 and our estimation of T^g. To be a cluster member 
in this work an object has to satisfy all four of these steps. 



3.1. First Candidate Selection: CMDs 



Candidates were first selected from our CMDs by keeping all objects which are no more 
than 0.15 mags redder or bluer than the isochrone in all CMDs (this number accommodates 
errors in the magnitudes and uncertainties in the model isochrone), plus errors from age 
estimation and distance to IC 2391 reflected on the isochrone. We additionally include 
objects brighter than 0.753 mag the isochrone in order to include unresolved binaries. In 
Figure [5] we show two CMDs for field 01 where candidates were selected based on 815/20 vs. 
815/20-914/27 and Rc vs. Rc-J {top 2 panels). We also present two CMDs from the deep 
field 32 using the medium band 770/19, 856/14 and 914/27 {lower 2 panels). From a total 
of 20 008 114 objects detected, 174 511 are kept (99.2% are rejected). 

We also present in this figu re low mass cluster member candidates from previous work 



which we de tected in our survey (jPatten fc Pavlovskylll999l . lBarrado y Navascues et al.ll2004 



Doddll2004l and X- ray sources detected by XMM-Newton where some are also presented in 
Marino et al.ll2005l ). Some c andidates fro n a prev ious studies are simply not detected in our 
work. This is the case with iPlatais et al.l (120071 ) where the faintest candidates have ^~15, 
which corresponds to 0.6 Mq (close to the saturation limit of the radial and outward fields 
at ~O.9M0 and at the saturation limit of the central deep fields, also at ~O.6M0). Also, no 
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objects in our sample match the 34 members studied by ISiegler et al.l (120071 ) because either 
our images are too deep so bright stars saturate (e.g. HD74275, HD74374 and VXR22a, 
which saturate in the short exposures), or the objects are not in our fields (e.g. VXR06, 
SHJMIO and PP07, which are between fields 32 and 37). Also, no objects match within 



4 arcsec between our objects and the 17 cluster candidates from iRoUeston fc Byrnd (119971 ) 
for similar reasons: the bright stars saturate (e.g. object ID 162, 311 and 362, which are 
birghter then our saturation limit in Rc), or the objects are not in our fields (e.g. object ID 
729 and 955, which are in the central part of the cluster between fields 20 and 27). 



Finally, we point out that the survey of IC 2391 based on proper motion done by iDodd 



(I2OO4J ) only covers an area of 1° diameter in the central part of the cluster. 



3.2. Second Candidate Selection: Colour-Colour Diagrams 

The second stage of candidate selection was achieved by taking all objects within 0.15 
mag of the isochrone of the NextGen model in selected colour-colour diagrams. In Figure 
O we present two colour-colour diagrams where only the objects from the first selection are 
plotted. Considering that many colour-colour diagrams are possible (we have 4-5 filters), 
the variation of colour as a function of T^s was used to reject the use of colours for which 
the NextGen model shows small variation in the M and L dwarf regime (this is illustrated 
in Figure [3 with the 815/20-914/27 colours). 

Because one source of contamination are backgro und red giants, we show theoretical 



colours for such objects using the atmosphere models of iHauschildt et al.l (jl999bl ). assuming 
that all objects have a mass of 5M0, 0.5 < log g < 2.5 and 2000 K < T^g < 6000 K. We can 
see that Rc-J vs. -Rc^815/20 is not best suited for selecting candidates since the isochrone is 
overlapped by red giant contaminants. However, in 815/20-J vs. 914/27- J, we see a clear 
distinction between the isochrone and the red giant contaminant in the brown dwarf regime 
(by more than 0.2 mag). This procedure definitely helps to remove red giant contaminants, 
and is further discussed in subsection 16. 3[ From a total of 174 511 objects, 33 794 are kept 
(80.6% are rejected). 



3.3. Rejection of Contaminants Based on Proper Motion 

Although the RMS error of our astrometry is 0.15-0.20 arcsec (WFI) and 0.3-0.4 arcsec 
(CPAPIR), we nonetheless estimated proper motions per year using the motion between the 
1999/2000 WFI data and the 2007 CPAPIR data in an attempt to reject objects which de- 
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viate significantly from tlie mean cluster proper motion in tlie literature. The typical error 
on our proper motion measurement is 24 milliarcsec per year (mas yr~^). The values of 
{fiaCos6,fis) for IC 2391 in mas yr~^ from the literature are (-25.04±1.53,+23.19±1.23), (- 
25.05±0.34 +22.65±0.28), (-24.64±1.13 , +23.25±1.23) and (-25.0 6±0. 25,+22.73±0.22) from 



Doddlfl2004[ ). lLoktin fc Beshenovl (120031 ) . ISanner fc Geffertl (120011 ) and lRobichon et all (119991 ) 
respectively. For our selection procedure we use the average of these, (-25.0±2.0,+23.0±1.7). 

We first investigated whether the cluster itself could be identified in the proper motion 
plane. To do this, we retained only those objects detected from our observation runs with 
WFI (1999, 2000 and 2007) and CPAPIR (2007) which have a match within 1 arcsec. We 
then examined the distribution in the {fiaCos6:fis) plane for any feature at (-25.0:+23.0). 
However, we see no clump in the distribution of the proper motion (Figure [H]). Considering 
the large errors and the absence of any structure at the expected location, we decided not to 
perform any selection using the proper motion of IC 2391. However, astrometry is used to 
remove all objects with a proper motion higher than 72 mas yr~^ (30") away from the cluster 
proper motion. 



3.4. Rejection of Objects Based on Observed Magnitude vs. Predicted 

Magnitude Discrepancy 

As indicated in §2.41 our determination of T^s is based on the energy distribution of each 
object and is independent of distance. The membership status is determined by comparing 
the observed magnitude of a given object in a band with the magnitude predicted based on 
its derived T^s and IC 2391 's isochrone. (The premise is that the predicted magnitude of a 
background contaminant would be lower - brigher - than its observed magnitude and higher 
for a foreground contaminant.) In order to avoid removing unresolved binaries that are real 
members of the cluster, we keep all objects with a computed magnitude of up to 0.753 mag 
brigher than the observed magnitude. In this procedure, we are also taking into acount 
photometric errors and uncertainties in the age and distance determinations of IC 2391. 
This is represented in Figure [91 Combined with the rejection of contaminants based on 
proper motion, this selection step reject 89.2% of the 33 794 candidates obtained from the 
CMD and colour- colour diagrams. 
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RESULTS OF THE SURVEY 



The final selection gives us 954 photometric candidates for outward fields (namely fields 
43, 46, 47, 48 and 49), 499 photometric candidates for the four deep fields (15, 20, 27 and 
32, with filters R^, 770/19, 815/20, 856/14 and 914/27) and 1 734 for all other radial fields 
(observed with filters Rc, 815/20, 914/27 and J). (We present in §5.31 a discussion of the 
contamination of the radial fields.) All our photometric candidates are presented in Table 
[21 Objects are given the notation IC 2391-WFI-ZZ-YYY where ZZ is the field number and 
YYY a serial identification number (ID). Only the first 10 rows of the tables are shown, 
the remainder available online. We also comp are in Table B] all objects in our sa mple which 
are al so confirmed as cluster members from iBarrado y Navascues et al.l (120041 ) and iDodd 
(120041 ) and was detected by the X-ray Multi-Mirror Mission ( XMM-Newton) . We see a 



good agreement between Tes from our photometric data and from lBarrado y Navascues et al. 
(120041 ). where only the colour {R — I)c was used to compute Tes- 



Not all candidate members from previous studies, detected in our survey, are mem- 
bers of IC 2391 based on our photometric selection. As pointed above, a cluster member 
presented in this work is an object that satisfies all four steps of our selection procedure. 
For instan ce, among the two o b jects detected in our survey which are also cluster can- 
didates by [Fatten fc PavlovskyI (Il999l ). one is recovered by our selection (object number 
8). From 10 objects classifi e d as c andidate members from spectroscopy and photometry by 



Barrado y Navascues et al.l (120041 ). 5 are recovered in our selection: objects CTIO-038 and 



-091 fail the colour-colour diagrams test and objects CTIO-041, -049 and -091 fail the pre- 
dicted magnitude v s. observed magnitude test . One possible source of disagreement could 
be the use of Ha by lBarrado y Navascues et al.l (120041 ) as a membership criteria (see §6.1l for 
further discussion of this issue). 



A total of 53 objects classified as cluster members by lDoddl (120041 ) based on proper mo- 
tion and photometry were detected in our survey. However, one of Dodd objects is recovered 
in our survey (object number 155, also identi fied as CTIO-152). This is also th e only matchs, 
within 4 arcsec, with the cluster members of iBarrado y Navascues et al.l (120041 ). Considering 
the size of the window used for their proper motion selection (in milliarcsec, -28>/iQ,cos5>- 
20 and +20</i5<+28) and the order of the error on the known proper motion of IC 2391 
(<2 miliarcsec in /i^cos^ and fis, see §3.31 below), one could suspect some contamination by 
fiel d star s . Thi s can be confirmed by the large scatter in the CMDs of IC 2391 presented 
by iDoddl (120041 ) below i?~12 and i^~10 in Figure 3 and 4 r espectively. Anothe r survey on 
IC 2391 based on proper motion was performed recently by jPlatais et al.l (120071 ). However, 
as discussed in §3.11 no objects from this work were detected in our survey. 



Only two objects listed Table 2 from 



Marino et al 



(12005! ) (IC 2391 members observed 
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with XMM- Newton) were also detected by our survey (source number 86 and CTIO-130), 
but neither objects is recovered by our photome t ric se lection. The first one is a source 
that overlap^ with VXR53 from Patten fc Simon ( 1996 ) and was identified as a suspected 
cluster member, and also overlaps with CTIO-126 from iBarrado y Navascues et al.l (120011 ) 
and was clas sified as a cluster member (howe ver, there was no spectroscopic follow-up of 
CTIO-126 by lBarrado y Navascues et al.ll2004l ). This object is not recovered in our selection 
because this candi date fail the observe d magnitude vs. predicted magnitude test. Another 
object presented by lMarino et al.l (120051 ) (observed only by the MOS cameras onboard XMM- 
Newton, but not by the pn camera) is CTIO-130, but they noted that this star has J and 
[J-K) values incompatible with the IC 2391 main sequence. 



4.1. Effect of Background Contamination on Candidate Selection 

In comparing the CMDs for different fields, we discovered something peculiar (Figure 
[TUl top 3 panels). We see a shift in the colours of the bulk of the (field) stars from field 
to field, something we also observe in other colours. The comparison of the amplitude of 
this shift (for a given magnitude interval) with observational parameters such as nights, 
airmass, seeing and lOcr detection limit shows no correlation and there is no other indication 
of reduction problems. We did, however, find a correlation of the colour shift with the 
Galactic longitude b. However, in order to verify that these shifts where real, we obtained 
DENIS photometry {Deep Near Infrared Survey of the Southern Sky) in / and J band for 
the same fields presented in Figure [TOl which are field 01, 09 and 40. We can see that the 
shift in the colours of the bulk of the (field) stars is also observed in the DENIS data (Figure 
[TUl lower 3 panels). 

Although reddening is negligible for objects in IC 2391, this is not the case for back- 
ground objects, and these constitute most of the stars in our sample. D ue to the high 



variation of the background extinction in this direction of the Galactic disk (jSchlegel et al. 



19981 ) - the cluster is centered at /=270.4 b=-Q.9 - some variation in the CMD locus could 
be extinction-induced variations in the background stars. In Figure [11] (left), we plot the 
reddening E{B — V) in our fields against the median of the colour 815/20-914/27 (in a 
bin of magnitude of 15 < 815/20 < 16) for all our fields. The colours vary by as much as 
0.25 mag. To better illustrate the spatial variation of the background extinction. Figure [TT] 



^The EPIC cameras on XMM-Newton have an angular resolution of 6 arcsec. Two of the cameras are 
MOS (Metal Oxide Semi-conductor) CCD arrays (referred to as the MOS cameras) and on camera at the 
focus of this telescope uses pn CCDs (referred to as the pn camera). 
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{right) show s the position o f the fields of our survey overplotted with the E{B-V) extinc- 
tion map of ISchlegel et aLl (119981 ) . This colo ur gradiant of the background stars has not 



been reported in previous surveys of IC 2391 (jPodd 



Patten fc Pavlovskv 1999[). It can be expected that 



2004; 



Barrado v Navascues et al 



Barrado y Navascues et al.l (l200ll ) and 



2001 



Patten &: Pavlovsky Jl999h didn't observed such shift in colour since their survey cover a 



smaller area (2.5 and 0.8 sq. deg. respectively) of the sky compared to our 10.9 sq. deg. 
coverage. 



4.2. Mass Function 

The mass function, ^ (logioM) , is generaly defined as the number of stars per cubic parsec 
(pc^) in the logarithmic mass interval logioM to logioM + dlogioM. Here, we do not compute 
the volume of IC 2391 so instead we present a MF using the total number of objects in each 0.1 
logioM bin per 1 000 arcmin^, starting at the mass bin logioM=-1.65 (~0.02Mq). The mass 
functions computed here are all system mass functions since we don't make any corrections 
for binaries. We analyse the radial variation of the MF using the fields with photometry with 
the filters Rc, 815/20, 914/27 and J (Figure [T2|) . Mass functions were computed over three 
regions: for fields between 0.5° to 1.5° of the cluster center (which corresponds to 1.3 pc and 
3.8 pc respectively); for the annulus from 1.5° to 2.1° (which corresponds to 5.4 pc); for fields 
outside of 2.1'^. We have also computed a MF for all fields within 2.1° of the cluster center 
to help radial variation analysis and present this as our estimation of the MF for IC 2391. 
Furthermore, we have measured the MF for the five outward fields and also for the four deep 
fields (Figure [T3|) . 

Since the radial fields were also observed with in addition to the filters used for 
both the radial and outward fields, we use this to estimate the (additional) contamination 
in outward fields relative to the radial fields for each mass bin. To do so, we performed a 
photometric selection for our radial fields using only the filters available in the outward fields 
(815/20, 915/27 and J). We compared the MF computed from this with the MF from the 
outward fields and obtained, for each mass bin, the number of object that would have been 
rejected if we would had an additional -Rc-band observation. (Here we make the assumption 
that the true MF should be the same in the radial and outwards fields.) 

In Figure [13] (left panel) we present the uncorrected MF of the outward fields and the 
corrected MF of the outward fields (right panel). It is not possible to perform such corrections 



•^For reference, the core and tidal radius of IC 2391 estimated by iPiskunov et al. 
7.38 pc respectively, which corresponds to 0.35° and 2.89° from the cluster center 



2007 



are 1.2 pc and 
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for the deep fields. 

Useful (an d simple) parametrizations of the mass function include the power law of 



Salpeterl (Il955l ) and a lognormal 



^(logio M) = • exp 



(logioM-logipMo: 
2(t2 



|2n 



(4) 



where fc=0.086, mo=O.22M0 and a=0.57 was derived for the Galactic field by IChabrier 



20031 ). Fitting the lognormal mass function to our data for all fields within 2.1° of the cluster 



center, we obtain fc=10.7±3.2, mo=0.13±0.03 Mq and (j=0.46±0.07. This is overplotted in 
Figure [131 

If we assume that the lognormal fit of the MF describe the behaviour of the population 
in the mass range 0.02-0.9 in IC 2391, the total number of object expected is 3 985 for 
a total mass of 679 Mq. (In §5.31 we will discuss again the total number of object and total 
mass, following an estimation of the contamination for each mass bin in the MF of the radial 
fields.) 



We present in Figure [14] the MF for all fields within 2.1° of the cluster center and from 

other open clusters with s imilar ages (NGC 2547, ~30Myr; IC 4665, 28"'"'''^ g.e) Myr based 

on Li depletion boundary, Manzi et al. 20081) . We also show on Figure [T^ th e MF of IC 239 1 
as determined in previous work (i.e. from lBarrado y Navascues et al.ll2004l and lDoddll2004j ). 
All were normalized to the Galactic field star MF at 0.3 Mq. 



5. ANALYSIS AND DISCUSSION OF THE MASS FUNCTIONS 

In the following subsection, we discuss the mass function derived from the deep fields 
and from outward fields only. The other fields are used to study the radial variation of the 
MF and are subject of further discussion in the following two subsections. We complete this 
section with a discussion of the contamination rate by non-cluster members. 



5.1. Mass Function of the outward fields and of the deep fields 

Considering the fact that only three bands were used for the outward fields, and thus 
fewer constraints imposed, we expect that the number of photometric candidates would 
be larger per unit area than the other fields. The MF (Figure [131 left panel) shows more 
low mass objects (compared to the MF of the radial field, for masses below ~O.15M0), a 
similarity with the radial MF from 0.13 to 0.3 Mq, and again more stellar objects in the mass 
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range 0.5 to 0.8 M©. The corrected MF (Figure [131 right panel) shows a better agreement 
with the radial MF. 

The MF from the deep fields (Figure [T^ agrees with the mass function of the radial 
fields within 2.1° from cluster center in the mass range 0.05 to 0.1 Mq and above 0.2 Mq. 
However, there is more substellar objects below 0.05 M©. 



The rise of the MF for objects below 0.05 M© was also observed in IC 2391 by lBarrado y Navascues et al 



( I2OO4J ). In their work, the MF was computed with objects that were selected as cluster mem- 
bers based on i?c; Ici J 1 H and K photometry. Since their NIR photometry was taken from 
2MASS, no data are available for objects fainter than /c^l9 (lOa detection limit of 2MASS 
is at J~15.8). As a result, their selection for fainter objects was based on and Jc photom- 
etry only. In our case, although J band photometry is available fo r the outward fields, no 



Rc photometry is available. Thus iBarrado y Navascues et al.l (120041 ) used a relatively short 



baseline {Rc-h) for their selection, as did we in our outward fields (815/20, 914/27 and J), 
both of which are considerably shorter than the baseline we used in the radial fields (i?c, 
815/20, 914/27 and J). This situation is also observed in the MF of the deep fields (where 
Rc band photometry is available, but no J band photometry). Only the fields observed with 
Rc and J as well (i.e. a longer baseline) show no significant rise of the MF below 0.05 M© 
(Figure [131 left panel). Since no red giants were observed in our spectroscopic follow-up 
of the two deep fields 15 and 20 (§6|), we conclude that this increase is an artefact due to 
contamination by M-dwarfs. In this low mass regime (for objects with mass <0.05Mq), a 
long spectral baseline (including, for instance, Rc and J) is needed to efficiently remove con- 
taminations, as it allows a better determination of the energy distribution. This is confirmed 
when we compare the corrected and uncorrected MFs of the outward fields (Figure [T3|) . 

The rise in the MF over 0.5-1.0 Mf^ is observed in the outward fields but not in the 



deep fields (Figure [13]). iJeffries et al.l mm present the MF of the open cluster NGC 2547 
and also noticed a rise in the 0.7-1.0 M0 interval (Figure This rise, also observed in 
the luminosity function as a large peak at 12</c<14.5, they attribute to contaminating 
background giants. This would be consistent in the fact that we see this in the radial fields 
(see Figure [T^ but not in the deep fields (Figure [T31). Indeed, as we discuss later in §S1 
no red giants were found in our spectroscopic follow-up, confirming that t he use of medium 



bands such as 770/19, 815/20, 856/14 and 914/27, and theoretical colours (IHauschildt et al. 



1999bl ) are effective in removing background red giants. However, from the MF of the radial 
fields (including the outwards fields), the medium filters 815/20 and 914/27 alone, combined 
with wide band Rc and/or J, are less efficient at removing background giants. Therefore, 
the rise in the MF of the outward fields over 0.5-1.0 M© is related to the filters used, but is 
not a baseline issue. 
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We also observe on Figure [13] that two mass bins at ~0.11 and ~0.18 Mq are significantly 
high compared to the other mass bins in the mass range 0.05-0.3 Mq. Considering that this 
is not observed in the MF of the radial field and in the corrected/uncorrected MF of the 
outward fields, we suspect that this is due to the selection procedure related to the deep 
fields so no conclusion should be made based on these two mass bins. 



5.2. Radial Variation of the Mass Function at the Stellar and at the 

Substellar Regimes 

From a first glance at Figure [T^ we see that the two mass functions within 2.1° (the 
points in the two left-hand panels) are somewhat similar. While there is some differences, 
these are not very significant compared to the difference between their common MF within 
2.1° (plotted as the histogram) and the MF in the outskirts of IC 2391 (beyond 2.1°) as 
shown in the right-hand panel. Indeed, the MF for 6'>2.1° shows a significant deficiency of 
stellar objects from 0.1 to 0.3 M© (logioM=-0.6) compared to the MF from the inner part 
of the cluster, whereas outside of this mass range there is no significant change with radius. 
Although we observe a number of objects in the highest mass bins at 0.5-0.7 M© (logioM=- 
0.3 to -0.15), we concluded (in §5.ip that this range of masses is subject to significant 
contamination by red giants, so no conclusion should be drawn from the radial variation of 
the MF in this mass interval. 

In Figure [15] we present the cumulative mass function for the same three regions of the 
cluster presented in Figure [12] We again see the relative absence of objects from 0.1 to 
0.3 M© (logioM=-0.6) for the inner radii. A Kolmogorov-Smirnov test performed on these 
distributions indicates that there is only a 1.1-10^^ % probability of getting such a difference 
under the null hypothesis that the population at 6'>2.1° is the same as that at 6'<2.1°, thus 
reinforcing the suggestion that the mass functions are significantly different. 

To help the analysis of the radial variation as a function of mass, we also present on 
Figure [TH] the radial profile of IC 2391 using the radial fields for four different mass bins: 
M< 0.072 M©, O.O72M0<M<O.15M0, 0.15 M© < M < 0.3 M© and O.3M0<M. (It is not 
surprising that a fit for the most massive stars is not possible since the core radius is at 
0.35° and we do not have any radial field s closer than 0.8°.) For the second and third radial 



profile, we fit a King profile (IKingll 19621 ). where the fit give us a maximal number density 



at the center of 27.7+17:3 and 26.11^:3 members per 1 000 arcmin and a full width at half 
maximum of 1.39°]']o;23 (or 3.5pc^|!,;5) and 1.45°^q;24 (or 3.7 pc't^l'^), respectively. This is 
purely for illustration purposes, as it is clear that this is not a good model for this mass 
profile. 
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In Figure [12] we do not see any significant radial variation of tlie MF in tlie substel- 
lar regime. Tlie radial profile of the substellar population in Figure [16] also indicates no 
significant radial variation. On the other hand, we do see reasonable evidence for a radial 
variation for masses above 0.072 Mq. From this radial profile and the mass functions already 
discussed, we can conclude that the spatial distribution of the BD population is uniform 
compared to the stellar population (from 0.072 to 0.3 M©), which is more clustered within 



Kumar fc Schmejal (120071 ) also found the stellar population to be more clustered than 
the substellar population in the clusters IC 348 and Trapezium. This would favour the 
ejectio n scenario for forrning B Ds if the BDs have a higher velocity dispersion than the 
stars (IKroupa fc Bouvierl l2003l ). because the higher velocity from ejection creates a more 
uniform spatial dis t ributi on for the BDs compared to the stars. However, the two clusters of 
Kumar fc Schmejal (120071 ) are both younger than 3 Myr, while IC 2391 has an age of 50 Myr, 
some ~3 times older than its crossing time {tcross =17 Myr). We can expect that if BDs 
have a higher velocity dispersion than stars, then in an older cluster most of the BDs with 
velocity dispersion great er than the escape velocity could have escaped the cluster already 
(iMoraux fc Clarkel boosh . 



The homogeneous distribution of the substellar objects compared to the more clus- 
tered stellar population could be instead a signature of mass segregation through dynamical 
evolution or of primordial origin. We have indicated previously that mass s egregation via dy- 
nami cal evolution could occur on a tim escale of order one relaxation time (iBonnell fc Davies 
19981 ) , or even less (lAUison et al.ll2009l ) . Considering that the cluster is only three times older 
than its crossing time, which may be insufficient for significant dynamical evolution, it is 
difficult to make an inference on the BD formation mechanism from the radial MF variation 
given the uncertainty about to what extend the cluster has evolved dynamically. If it could 
be demonstrated that the cluster is dynamically unevolved and that the BDs have a higher 
velocity dispersion than the stars, then our observations are consistent with BD formation 
by the ejection hypothesis. 



5.3. Contamination by non-members 



Possible point so urce contaminants oth er than field M-dwarfs include red giants and 
high redshift quasars (ICaballero et al.ll2008l ). However, as pointed out in §5.11 and in our 
spectroscopic follow-up (see below in §6.31) . the use of medium band filters and theoretical 
colours is efficient at removing potential background red giant contaminants. As for the high 
redshift quasars (for instance at z^6), their spectral energy distribution is similar to mid-T 



- 19 - 



dwarfs and moreover, they are rare (ICaballero et al.ll2008l ). Considering that our faintest 



targets are early L-dwarfs, the MF should not be affected by contamination by quasars. 

Here we present an estimation of the contamination in our photometric survey based on 
the radial fields. First, we used the radial profile of IC 2391 using the radial fields for four 
different mass bin (Figure [T6l) . We assumed that near the tidal radius (~2.89°), the number 
of objects per 1 000 arcmin^ should be zero. From this, we computed that we can expect a 
contamination of ~8.8 objects per 1 OOOarcmin^ for masses above 0.3 M©, ~8.8 objects per 
1 OOOarcmin^ in the mass range of 0.15 < M < 0.3 M©, ~16.7 objects per 1 OOOarcmin^ 
in the mass range of 0.072 Mq < M < 0.15 and ~9.8 objects per 1 OOOarcmin^ in the 
substellar regime. If we use the same assumptions about the lognormal fit of the MF, the 
new total number of object expected in IC 2391 is 1 954 for a total mass of 308 Mq. 



6. PRELIMINARY SPECTROSCOPIC FOLLOW-UP 

Here we present the results of a preliminary spectroscopic follow-up of some photometric 
candidates. As explained in the previous section, the main sources of contamination in our 
photometric selection are background red giants and field M-dwarfs. We have also shown 
in §4. II that, because of extinction, background contamination is non-uniform. We therefore 
need to refute or confirm membership status with optical spectra. For this task we used 
the fiber spectrograph HYDRA. It is not possible to cross fibres with this instrument, so 
we have not yet been able to observe all candidates in a given field. (It is our intention to 
eventually obtain spectra of all candidates.) The data reduction was described in §2.51 while 
the spectral type and luminosity class determination was presented in ^2.61 The spectro scopic 



Teff was obtained using the spectral type and the temperature scales of lLuhmaru (119991 ) while 
each mass was derived from Teg using our isochrone for IC 2391. We discuss membership 
determination based on optical spectra below in §6.2[ 



Among the spectra obtained, 17 had a signal-to-noise ratio (SNR) higher than 5. These 
are presented in Figure [T71 Table H] provides the derived parameters (spectral type, T^s 
and mass) and the SNR. Objects are given the same notation as the photometric candi- 
dates : IC 2391-WFI-ZZ-YYY where ZZ is the field number and YYY a serial identifi- 
cation number (ID). Table [5] giv es details of an object confirmed as cluster members by 



Barrado y Navascues et al.l (120041 ). including their SpT and T^s determination, for which we 



also have a spectra. 
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6.1. Contaminating Ha nebula emission 



We mentioned in §2.51 the presence of contamination at Ha in sky spectra. As tliese 
spectra are used for background subtraction (we have fiber spectra), there are potential 
difficulties in measuring the stellar Ha line. We now discuss this issue. 

In produ cing a high-resolu ti on at las of night-sky emission lines with the Keck echelle 
spectrograph, lOsterbrock et al.l (119961 ) observed an Ha emission line at high Galactic lat- 
itudes which they concluded was due to diffuse interstellar gas emission (the closest at- 
mospheric emission observed were two OH lines at 6553.617A and 6568.779 A). From the 
AAO/UKST SuperCOSMOS Ha Survey fSHS. IParker et al.l boO,'^ ). we have also noticed 
high variations of Ha emission at low Galactic latitude. We used the SHS to estimate the 
Ha emission at each position of our sky fibers (by taking a median of the fiux over a 200 
X 200 arcsec window). The frames are fiat-field corrected but not fiux calibrated, so we 
retain the unit of (photon) counts. In Figure [18] we plot this against the fiux (in counts) 
of the Ha emission line of our background spectra for field 20. While there is no strong 
evidence for a correlation, we nonetheless see a significant variation of the Ha emission. As 
this clearly prevents a reliable background subtraction of the Ha line, we choose not to draw 
an y conclusions on membership statu s based on this line. We must therefore question its use 
by iBarrado y Navascues et al.l (120041 ) for this purpose (who used the same instrument for 
the same cluster). For further observations of objects in direction of IC 2391 using fiber-fed 
spectrograph, w e recommend backgrou nd subtraction to be performed in a similar way as 
the one done by ICarpenter et al.l (Il997l ). where the same fibers for the science targets were 
also used for sky subtraction but shifted 6 arcsec away. 



6.2. Membership Determination 



We use the Li I line at 6708 A to help confirm substellar status of photometric candidates 
and to establish membership of IC 2391. Lithium can be observed in young, more massive 
stars with radi ative interiors bec ause of less efficient mixing than in fully convective low 
mass stars (e.g. iManzi et al.ll2008l ). Lithium may still even be present in the atmospheres of 
young, fully convect ive low mass stars, if they are young enough that not yet all the lithium 
has been "burned" (IManzi et al.ll2008[ ). Older, lower mass BDs (<O.O65M0) never achieve 
core tempertures h igh enough to burn lithium and so preserve their initial lithium content 
( jRebolo et al.l 1 19961 ). Here we assume that field stars (with M>O.O72M0) are too old to still 
retain lithium in their atmospheres. Hence we take the presence of the Lil line in candidates 
with M<0.072 Mq as an indicator of membership in IC 2391, as only cluster members fa inter 
than the Lithium depletion boundary (Jc =16.2 mag. iBarrado y Navascues et al.ll200ll ) are 



- 21 - 



young enough to retain Lithium. 

The sodium doublet at 8200 A is a gravity indicator, so is someti mes used to exclude 



field stars. Specifically, its equivalent width (EW) is sensitive to \ogg ( iMartm et al.lll996l ) 
and because BDs contract as they age, logg will increase. Because field late M-dwarfs (which 
have similar colours to M dwarf cluster members) will generally be much older and so more 
evolved, they will have larger EWs in this li ne (for a given chemical compositi on). We use 



the EW measurement of CTIO-046 from the iBarrado y Navascues et al.l (12004 ) survey as a 



lower limit on EW values for M-dwarfs to be non-members. (This object was defined as a 
non-member based on various criteria and had W^(NaI)=7.3±0.2 A). It can be argued that, 
since surface gravity changes with mass, the limit EW(Nal)=7.3A could also change with 
mass. Based on our isochrone of IC 2391, from 0.04 to 0.2 Mq (which is the mass range of 
our spectroscopic follow-up), the surface gravity will varies only from log 5^ = 4.65 to 4.72 
(51ogg = 0.07). 

The spectral resolution is sufficient to provide an estimate of the radial velocitj]^. As 
for the radial velocity criteria, we exclude candidates which differ significantly (±3a) frona 



a recent determination of the cluster's radial velocity (16±3km/s, iKharchenko et al.l 12005 
where a is the the error of the radial veolicty of IC 2391 added in quadrature with the 
error of our candidates). We didn't used radial velocity measurment for which errors exceed 
30km/s, which is ten time the error on the radial velocity of IC 2391. 

Finally, we use the SpT determination to obtain Tes and masses for each spectrum. In 
order to be confirmed as a (spectroscopic) cluster member, the spectroscopic T^s and must 
agree with the photometric T^s to within 200 K. 

In Table E] we again present all objects from Table HJ but with physical parameters 
and with membership status based on photometry and spectroscopy (i.e. which satisfy our 
spectroscopic criterium). We don't reject objects below which does not present a feature of 
Lil due to low SNR if the other criteria are satisfied (e.g. IC2391-WFI-15-005). 



6.3. Discussion of the spectral data 

Of the 17 photometric candidates observed with a SNR higher than 5, 9 are spectroscopic 
members of the cluster. We find no red giants in our spectral sample, which demonstrates 
that our choice of filters and selection procedure is efficent at minimizing this contamination. 



"^The radial velocity measurment was performed with the IRAF task xcsao. This task perform cross 
correlation against a spectrum with known radial velocity and makes the barycentric correction. 
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Since our spectroscopic follow-up covers only part of the mass range used for the mass 
function calculation in §5l it is not possible to compute a new MF with corrections applied 
at each mass bin. It is expected that the contamination rate would be different for other mass 



ra nge and fields (for different f i 



by iBarrado y Navascues et al- 



ter co mbinations). One of our spectral targets was observed 



mm-. CTIO-62, which has the label IC2391-WFI-20-067 in 



Barrado y Navascues et al.l (120041 ) on the status (cluster member) 



our survey. We agree with 
of this object. 

Another of our spectral targets, IC2391-WFI-20-001, shows Lithium even though this 
object has an inferred mass above the stellar/substellar boudary (spectroscopic mass of 
0.079 Mq and a photometric ma ss of 0.089 Mfr^). The lithium dep letion boundary has been 
estimated to be at Ic =16.2 mag (IBarrado y Navascues et al.ll200ll ). which corresponds to an 
effective temperature of about 3 000K based on our isochrones of IC 2391. As the cluster is 
not that old, the lithium depletion boundary lies well above the substellar boundary, so its 
not surprising to see a trace of Lithium in the spectrum of this object (which has a magnitude 
of Ic =16.681 mag). 

Although we obtained for IC2391-WFI-20-024 the same spectral type as for IC2391- 
WFI-20-029, we consider this object as a member of IC 2391, but to have a mass slightly 
above the substellar. Its photometry gives Teg = 2958 K and M = 0.081 M©, with a predicted 
magnitude of 815/20 = 16.445 mag. If this object were an unresolved binary, we would expect 
its observed magnitude to be brighter than its predicted one, which is not the case (observed 
magnitude of 815/20 = 16.520 mag). Also, it should be noted that this object was observed 
with a fibre with poor spectral response below 6800 A (which is why we don't show it in 
Figure [T71). Although this does not affect the PC3 index used for the SpT determination 
(which covers 7540 A to 7580 A, and 8230 A to 8270 A), it does influence the reduction pro- 
cess (including the throughput correction, illumination correction, extraction of the spectra 
and flux calibration). Therefore, we consider its SpT uncertainty to be larger (two subtypes 
rather than one). 



6.4. Discovery of new brown dwarf members of IC 2391 

Of the 17 spectral targets, we assign as brown dwarfs two in IC 2391, on the basis of 
spectroscopic confirmation, and having both photometric and spectroscopic masses below 
O.O72M0. These are new discoveries. These objects are IC2391-WFI-15-005 and IC2391- 
WFI-20-029. Table [7] lists their parameters. Figure [19] shows their spectra and Figure [20] 
contains the finding charts. We can see in Figure [T9] th at Ha is not visible in IC239 1 -WFI - 
20-029, so it would be designated as non-members by IBarrado y Navascues et al.l (12004] ). 
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Considering that the MF from our radial fields is similar to that of the deep fields in the 
mass range of these two new objects (from 0.045 to O.OTMq), then if we used the same 
selection method, then statistically we would expect to find two brown dwarfs in the same 
mass range in the other two deep fields, and about 31 in all of the radial fields. 



7. CONCLUSIONS 

We have performed a multi-band photometric survey over 10.9 square degrees of the 
open cluster IC 2391, and completed a preliminary spectroscopic follow-up of brown dwarfs 
and very low mass stars candidates from two of the WFI fields. Our objective was to study 
the mass function of this cluster, and in particular its radial dependence. We observed a 
radial variation in the MF from 0.072 to 0.3 M©, but we do not observe a significant radial 
variation in the mass function in the substellar regime. This comparative lack of radial 
variation of the substellar mass function is in favour of the ejection scenario for forming 
brown dwarfs, but considering that IC 2391 is ~3 times older than its crossing time, we 
might expect that most of the brown dwarfs with velocity dispersion greater than the escape 
velocity could have already escape the cluster. On the other hand, the rather homogeneous 
distribution of the substellar objects and the clustered distribution of stellar objects within 
9r^2° could be a signature that mass segregation via dynamical evolution has occurred in 
IC 2391, or that this mass segregation is of primordial nature. We have concluded that if this 
cluster is dynamically unevolved and if the brown dwarfs have a higher velocity dispersion 
than the stars, then our observations are consistent with brown dwarf formation by the 
ejection hypothesis. 

In addition to the radial study, we derived a mass function from four central deeper 
fields as well as from five fields near the edge of the cluster observed with only three filters 
(the outward fields). In both cases we see an apparent rise in the number of objects below 
0.05 Mq (logioM=-1.3), but we concluded that this is an artefact of residua l contamination 



by field M dwarfs. This was also seen by iBarrado y Navascues et al.l (120041 ). The fact that 



we don't see this rise in the radial fields is because they were observed with both the J 
and Rc filters in addition to the medium band filters. This longer spectral baseline permits 
a better determination of the energy distributions and thus helps the rejection of objects 
(in particular field M dwarfs) based on observed magnitude vs. predicted magnitude from 
models. 



Another apparent rise in the MF over the 0.5-1.0 M© interval (also observed by I Jeffries et al 



2004J for NGC 2547) is due to background giants. Red giant contamination may be reduced 



by using medium bands such as 770/19, 815/20, 856/14 and 914/27, and theoretical colours 



of red giants (IHauschildt et al.lll999bl ). Our spectroscopic follow-up has confirmed that se- 
lection based on these filters resulted in no red giant contaminants among our sample of 
spectra. 



We see some variation in the colours of the main (field star) locus which we attribute to 
variable extinction affecting the background stars. This underlines the need for spectroscopic 
observations in this cluster in order to confirm membership and/or brown dwarf status in 
individual cases. 

We have performed a preliminary spectroscopic follow-up of photometric cadidates in 
two of our deep fields (0.5 sq. degrees). Of 17 photometric candidates, we confirm 9 objects 
(i.e. half) as true cluster members. Of these, two are new brown dwarf members of IC 2391 
(in the sense that they fufill our spectroscopic and photometric criteria). Using our derived 
mass functions for the deep and radial fields, we expect there to be two more brown dwarfs 
in the mass range 0.045 to 0.07 in the other deep fields and up to 31 in all the other 
radial fields in the same mass range. 

Finally, we find that the Ha line cannot be used as a membership criterion from fiber 
spectroscopy at low spectral resolution (spectral dispersion of 1.14Aper pixel) because of 
spatially variable diffuse Ha emission. This prevents reliable sky subtraction around this 
line when using a fiber spectrograph with fibers assigned for sky subtraction. 
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Fig. 1. — Area covered by our survey of IC 2391. The four thick squares are the deep fields, 
the dotted squares are the outward fields and the others are the radial fields. The circles 
have radii of 1°, 2° and 3° from the cluster center. 
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Fig. 2. — Transmission curve of the filters used in our survey compared to the synthetic 
spectrum of a brown dwarf with = 2300 K, log g = 4.5 and solar metallicity (NextGen 
model). 
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Fig. 3. — Estimation of the completeness limit for the radial part of our survey using the 
J band (left) and of the deep part using 914/27 (right). The tick lines give best linear fit 
before the turn off; the vertical tick dotted line is the 10 a detection limit and the vertical 
thin line is the magnitude for which saturation start to occur in the short exposures. 
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Fig. 4. — Spectrum used for sky subtraction of our spectroscopic data. Note the Ha (nebula) 
emission line with equivalent width of iy(HQ;) =48 A. 
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Fig. 5. — Top. Two CMDs from the radial field 01. As dotted lines we show the isochrone 
computed from an evolutionary model with a grainless atmosphere (NextGen model, the 
masses for each are shown in the right panel). The thin dashed li ne is the lOcr detec- 
tion limit. We also show candidate low mass cluster members from Patten fc Pavlovsky 



(I1999I ) (gtars) ■ iBarrado y Navascues et al.l (12004J ) (triangles) . \Do(iai (12004J ) (squares) and X- 
ray sources detected by XMM-Newton (circles), which we detected in our survey. Some of 
these objects are not present in the left panel since the deep fields, where most of these 
objects are detected, lack J-band photometry. Bottom. Two CMDs from the deep field 32 
with medium bands. Isochrones, lOa detection limit and cluster members from previous 
studies are the same as for the top two panels. In each panel, the thin solid lines represent 
the selection curve and the errorbars present the photometric errors. 
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Fig. 6. — Two colour-colour diagrams of objects that are candidates based on o u r first 
selection (f rom fie l d 01) . Isoclirones and objects from iBarrado y Navascues et al.l (120041 ) 
{triangles) , iDoddl (120041 ) (squares) and XMM-Newton (circle) are as shown in Figure O 
The thin lines represent the colour of possible background red giant contaminants. The 
colour-colour diagram on the left is therefore not suited for further candidate selection since 
the isochrone is spanned by the colour grid of the red giants (which is not the case for the 
colour-colour diagram on the right). In each panel, the thin solid lines next to the isochrone 
represent the selection curve. 
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Fig. 7. — Colour sensitivity to effective temperature for the medium band filter 914/27. The 
vertical lines represents are at the approximate effective temperatures for spectral classes 
M5V (dash-dotted fine) and K5V (long-dash line). We can see that there is no variation at 
all for 815/20-914/27 compared to other colours at the LOV and M5V regime. 
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Fig. 8. — Left. Proper- motion diagram from our survey (in milliarcsec per year). IC 2391 
is at (-25.0, +23.0), which is taken from the literature ( §3.3p . The typicaly error bar of 
individual objects in our survey is shown. For clarity, only one object out of five is shown. 
Right. Contour plot of the same data. 
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Fig. 9. — Difference l:)etween the observed 8 f 5/20 magnitude and tliat computed from tlie 
derived mass and T^s, as a function of Tgg. Tlie four vertical lines are at the positions of 
LO, M5, K5 and G5 dwarfs (left to right). The dotted line (at —0.753) represents the error 
due to the possible presence of unresolved binaries at, the dashed-dotted lines represents the 
error on the magnitude determination and the long dashed lines represents the uncertainties 
on the age and distance of IC 2391. (The short-dashed line just traces zero). 
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Fig. 10. — CMDs of 3 fields from our survey {top, from left to right, fields 40, 01 and 09) 
and 3 CMDs from the same fields using DENIS data {bottom, I versus I- J). The NextGen 
isochrones is also shown. We clearly see a colour shift of the (field star) locus between these 
fields in our data, as well in the DENIS data. For all panels, the arrows represents the 
reddening vectors based on E{B-V) towards each fields. 
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Fig. 11. — Left. E{B — V) towards all our fields (from the lSchlegel et al.lll998l extinction map) 
plotted against the median 815/20-914/27 stellar colour in those fiel ds. Right. Pos i tion o f 
the fields of our survey overplotted by the E{B-V) extinction map of ISchlegel et al.l (1l998l ). 
The contour separation is 0.2 mag. 
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Fig. 12. — MF based on photometry for all radial fields. The lOcr detection limit is shown 
as a vertical dashed hne. Dots in each panel represent the MF of (left) fields within 1.5° 
of the cluster center, (center) fields within the annulus from 1.5° to 2.1° and (right) the 
MF of fields outside of 2.1°. Error bars are Poissonian arising from the number of objects 
observed in each bin. The histogram is the MF for all fields within 2.1° of the cluster center. 
The vertical thin dotted line is the mass for which saturation start to occur in the short 
exposures. (The total area covered for each panel, from left to right, is 6 637, 9 539 and 
5 609 arcmin^.) Just for reference, the ordinate value of 1.13 for the the bin at logioM=-0.85 
(0.14 Mq) - the histogram peak - corresponds to 109 objects. 
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Fig. 13. — Left. Filled dots represent the MF based on the four deep fields (observed with 
the wide bands Rc and the medium band 770/19, 815/20, 856/14 and 914/27) and open 
dots represent the MF based on the outward fields (observed with the wide band J and the 
medium bands 815/20 and 914/27). Also, we present the lOcx detection limit, the MF of 
all fields observed in 815/20, 914/27 and J within 2.1° from the cluster center and its 
log normal fit. The vertical thin dotted and thin dashed fine fines are the mass for which 
saturation start to occur in the short exposures for outward and deep field respectively. 
Right. Same as the left panel, but the outward fields MF presented is the corrected MF. 
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Fig. 14. — MF of IC 2391 from our present work (em yty dots, froni all fie lds within 2.1° from 



the cluster cente r), for IC 4665 (filled triangles) from lde Wit et al.l (120061 ) and for ^ 
(filled dots) from i Jeffries et al.l (120041 ). We also show the galactic f ield stars MF from 



(2003) as a dashed line. We also present the MF of IC 2391 from lBarrado y Navascues et al 



GC 2547 



Chabriei 



( I2OO4J ) (empty triangles) and from from iDoddl (120041 ) (empty squares). 
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Fig. 15. — Cumulative number of cluster members within 1.5° (dash-dotted line), within the 
annulus from 1.5° to 2.1° (dotted line) and outside of 2.1° (tick line). The lOo" detection 
limit is shown as an horizontal dash line. The vertical thin dotted line is the mass for which 
saturation start to occur. 
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Fig. 16. — Radial profile of IC 2391 for four different mass bins as indicated in each panel. 
A King profile fit JKing|ll962[ l was used for the radial profile for 0.072 Mq <, M < 0.15 Mq 
and 0.15 Mq <, M < 0.3 Mq (solid line). An estimation of the background contamination for 
each mass bin is given by the horizontal dotted line. In the panel of 0.15 M© <, M < 0.3 Mq, 
the four data points at N = indicate that no object was detected in that mass range for 
these fields (fields 03, 04, 31, 40 and 41, where 03 and 41 are at a similar distance of ~2.94° 
from cluster center). 
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Fig. 17. — Spectra used in our analysis. 
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Fig. 18. — Left. Ha emission (median counts) in fields of size 200 x 200 arcsec from the 
SuperComsos Survey at the location of our sky fibers, plotted against the flux (in counts) 
of the Ha emission hne we measured in our sky fibers in the HYDRA pointing of field 20. 
Right. Ha observations from SuperCOSMOS towards field 20 (images of 35 x 35 arcmin) . 
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Fig. 19. — Spectra of the eight newly discovered brown dwarf members of the IC 2391 cluster 
found in our survey. We also present, for each spectra, a new close-up on the Li and the Nal 
doublet. 
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Fig. 20. — Finding charts of the two new brown dwarf members of IC 2391 (from 815/20 
images). The panels are 3.5x3.5 arcmin with North up and East to the left. 



Table 1. Description of observations. 



Field 




RA 


DEC 


Distance (°) 


Region name 


Rc 


770/19 


815/20 


856/14 


914/27 


J 


01 


8 


24:38.8 


-51:18:16.5 


2.966 


radial 


1500/23.5 




1800/20.8 




600/19.7 


1820/17.8 


03 


8 


28:10.8 


-50:46:50.0 


2.945 


radial 


1500/23.8 




1800/20.5 




600/19.7 


1820/17.7 


04 


8 


27:56.8 


-52:07:47.9 


2.084 


radial 


1500/23.9 




1800/20.5 




600/19.5 


1820/17.5 


05 


8 


29:06.5 


-52:26:14.6 


1.793 


radial 


1500/23.7 




1800/21.1 




600/19.8 


1820/17.5 


06 


8 


30:30.6 


-51:38:59.1 


2.051 


radial 


1500/23.8 




1800/21.1 




600/19.7 


1820/18.0 


08 


8 


32:01.8 


-52:15:26.1 


1.481 


radial 


1500/23.9 




1800/21.1 




600/19.9 


1820/17.9 


09 


8 


33:15.7 


-50:39:24.6 


2.640 


radial 


1500/23.9 




1800/20.8 




600/19.8 


1820/18.0 


10 


8 


33:20.3 


-51:50:14.6 


1.616 


radial 


1500/24.0 




1800/21.0 




600/20.0 


1820/18.0 


11 


8 


34:06.5 


-51:24:35.9 


1.904 


radial 


1500/23.5 




1800/20.7 




600/19.8 


1820/17.9 


12 


8 


34:02.3 


-52:45:43.2 


0.978 


radial 


1500/23.7 




1800/21.1 




600/20.0 


1820/17.9 


14 


8 


36:17.4 


-50:38:44.9 


2.498 


radial 


1500/23.4 




1800/20.8 




600/19.7 


1820/17.9 


15 


8 


38:31.6 


-53:35:29.4 


0.757 


deep 


3900/22.7 


3900/20.9 


3000/20.7 


1500/19.7 


3000/20.3 




17 


8 


37:37.6 


-51:34:05.6 


1.552 


radial 


1500/23.6 




1800/20.9 




600/19.7 


1820/17.5 


18 


8 


38:11.5 


-52:01:44.1 


1.085 


radial 


1500/22.7 




1800/20.9 




600/19.5 


1820/17.6 


19 


8 


38:17.7 


-50:58:08.4 


2.121 


radial 


1500/24.0 




1800/20.6 




600/19.6 


1820/17.7 


20 


8 


38:45.1 


-52:35:58.0 


0.519 


deep 


3900/22.4 


3900/20.8 


8400/21.3 


15678/21.1 


3000/20.3 


~ 


21 


8 


41:22.5 


-52:14:04.7 


0.854 


radial 


1500/23.8 




1800/21.2 




600/19.9 


1820/17.7 


22 


8 


40:50.6 


-51:31:11.6 


1.553 


radial 


1500/23.7 




1800/21.1 




600/20.0 


1820/17.8 


24 


8 


39:59.7 


-54:14:24.0 


1.170 


radial 


1500/23.8 




1800/20.5 




600/19.7 


1820/17.9 


26 


8 


40:16.2 


-55:16:12.0 


2.200 


radial 


1500/23.7 




1800/20.6 




600/19.7 


1820/17.7 


27 


8 


41:01.5 


-53:50:51.7 


0.789 


deep 


3900/22.5 


3900/20.7 


9300/21.5 


7800/20.4 


4800/20.7 




28 


8 


41:46.0 


-54:46:27.8 


1.720 


radial 


1500/23.8 




1800/21.1 




600/19.8 


1820/17.5 


31 


8 


44:09.0 


-55:27:58.0 


2.464 


radial 


1500/23.6 




1800/20.5 




600/19.5 


1820/17.5 


32 


8 


44:10.1 


-52:39:49.0 


0.721 


deep 


3900/22.3 


3900/20.9 


8450/21.4 


10500/20.7 


4800/20.5 




35 


8 


44:39.7 


-54:21:53.7 


1.453 


radial 


1500/23.5 




1800/21.0 




600/19.9 


1820/17.2 


37 


8 


45:54.4 


-53:25:53.5 


0.927 


radial 


1500/23.5 




1800/21.1 




600/19.7 


1820/17.5 


38 


8 


47:00.8 


-53:55:12.1 


1.323 


radial 


1500/22.5 




1800/20.9 




600/19.6 


1820/17.7 


40 


8 


48:02.4 


-55:09:16.4 


2.380 


radial 


1500/22.8 




1800/20.6 




600/19.6 


1820/17.7 


41 


8 


48:09.8 


-55:46:27.3 


2.942 


radial 


1500/22.7 




1800/20.4 




600/19.4 


1820/17.7 


42 


8 


49:10.6 


-54:35:58.2 


2.023 


radial 


1500/22.7 




1800/21.1 




600/19.9 


1820/17.5 


43 


8 


50:15.8 


-53:23:37.1 


1.540 


outward 






1800/20.5 




600/19.3 


1820/17.5 


46 


8 


53:03.1 


-54:23:52.2 


2.318 


outward 






1800/20.5 




600/19.6 


1820/17.3 


47 


8 


53:26.1 


-53:52:29.4 


2.127 


outward 






1800/20.9 




300/19.3 


1820/17.3 


48 


8 


54:39.8 


-53:29:44.4 


2.203 


outward 






1800/20.5 




600/19.7 


1820/18.2 


49 


8 


56:47.5 


-54:17:22.5 


2.742 


outward 






1800/20.7 




600/19.6 


1820/17.8 



Table 1 — Continued 



Field RA DEC Distance (°) Region name Rc 770/19 815/20 856/14 914/27 J 



Note. — System notation is exposure time in seconds / 10a detection limit while — indicate that no 
observations was performed for that field in that filter. Distance give the distance of the field from cluster 
center (in degree). The lOcr detection limit in i?c is smaller than for the deep fields, although the exposure 
time was higher. Rc observations of the deep field were done in January 2000, while the observations of 
the radial fields were done in April 2007. Condensation problems and non-photometric nights was reported 
for the January 2000 observation run. 



Table 2. All photometric candidates of our survey 



Field 


ID 




RA 


DEC 




770/19 


815/20 


856/14 


914/27 


J 


M 


Tcs 


[815/20J 


01 


001 


8 


26:16.055 


-51:02:52.32 


19.645 




16.706 




16.292 


15.023 


0.057 


2768 


17.220 


01 


002 


8 


26:14.211 


-51:02:47.64 


18.426 




16.200 




15.920 


15.145 


0.104 


3072 


15.943 


01 


003 


8 


25:28.826 


-51:13:56.71 


18.866 




16.383 




15.998 


14.909 


0.081 


2957 


16.450 


01 


004 


8 


25:47.698 


-51:06:25.26 


18.848 




16.332 




15.953 


14.969 


0.078 


2935 


16.537 


01 


005 


8 


25:26.432 


-51:03:59.27 


18.359 




16.077 




15.771 


14.962 


0.096 


3043 


16.080 


01 


006 


8 


25:01.148 


-51:03:37.82 


19.574 




16.953 




16.582 


15.446 


0.072 


2890 


16.710 


01 


007 


8 


24:10.028 


-51:07:26.14 


19.283 




16.735 




16.341 


15.427 


0.076 


2925 


16.576 


01 


008 


8 


24:33.508 


-51:05:35.22 


18.764 




16.437 




16.094 


15.130 


0.093 


3029 


16.146 


01 


009 


8 


23:19.415 


-51:32:35.82 


19.900 




17.106 




16.655 


15.618 


0.065 


2838 


16.930 


01 


010 


8 


23:29.455 


-51:32:00.48 


18.373 




16.076 




15.762 


15.014 


0.094 


3032 


16.134 



Note. — Table [2] is published in its entirety in the electronic edition of the Astrophysical Journal. A portion is shown here 
for guidance regarding its form and content. The error on the determination of masses and effective temperature are the follow- 
ing : ATcff = 140K and AM = 0.1Mq for stars (M>O.2M0), ATefi =230K and AM = 0.05Mq for VLMS (0.072 q <M<0.2Mq), 
AToff = 420 K and AM = 0.02 M0 for BDs (M < 0.072 Mq). The magnitude [815/20] is the predicted magnitude based on photometric 
determination of Tgff and mass. 



Table 3. Objects from IBarrado y Navascues et al.l (120041 ). iDoddl (120041 ) and detected by XMM- Newton which are 

photometric candidates in our sample 



Field ID RA DEC 815/20 M T^s [815/20] NAME Ic {R-I)c T^s 



18 


006 


8 


38:47.074 


-52:14:56.16 


17.076 


0.053 


2723 


17.396 


CTIO-061 


17.309 


2.141 


2801 


20 


028 


8 


38:47.282 


-52:44:32.61 


16.662 


0.077 


2927 


16.567 


CTIO-062 


16.765 


2.000 


2937 


27 


002 


8 


40:09.537 


-53:37:49.81 


16.153 


0.095 


3036 


16.115 


CTIO-077 


16.308 


1.929 


2960 


32 


120 


8 


44:02.109 


-52:44:10.73 


17.050 


0.065 


2842 


16.911 


CTIO-160 


17.151 


2.090 


2806 


32 


231 


8 


43:38.421 


-52:50:55.15 


14.810 


0.206 


3310 


14.707 


8 


14.890 


1.790 




32 


295 


8 


43:38.422 


-52:50:55.13 


14.848 


0.206 


3310 


14.707 


CTIO-152 


14.891 


1.781 


3053 


32 


295 


8 


43:38.422 


-52:50:55.13 


14.848 


0.206 


3310 


14.707 


155 


14.530 


2.260 




32 


325 


8 


46:15.404 


-52:49:37.61 


15.305 


0.170 


3249 


15.046 


2XMM J084615.3-524937 








32 


340 


8 


46:04.238 


-52:45:18.99 


15.921 


0.122 


3134 


15.640 


2XMM J084604.3-524518 








37 


024 


8 


47:07.572 


-53:09:45.32 


15.377 


0.159 


3228 


15.164 


2XMM J084706.2-530944 









Not e. — Object 8 is from lPatten fc Pavlovskvl l|l999h . CTIO objects are from lBarrado v Navascues et al" ] ||2004h . object 155 is from lDoddl 
l l2004h while the 2XMM objects are from XMM-Newton. The values of 7c, [R — I)c and T^ff given in the last three columns are taken from ^ 
each previous studies. I 



Table 4. Stellar parameters from the objects observed in our spectroscopic follow-up. 



Field 


ID 




RA 


DEC 


SNR 


SpT 


-'eft 


M 


81 5/20 


-'eft ^,ptlut^ 




[81 5/201 


15 


003 


8 


36:40.363 


-53:21:30.00 


6.1 


M2.0 


3510 


0.413 


16.201 


3020 


0.091 


16.191 


15 


004 


8 


36:04.075 


-53:29:25.00 


20.7 


M3.5 


3265 


0.173 


16.306 


3010 


0.089 


16.232 


15 


005 


8 


36:18.241 


-53:25:57.60 


8.3 


M7.5 


2660 


0.049 


17.803 


2654 


0.048 


17.672 


15 


010 


8 


35:21.606 


-53:42:05.12 


8.0 


M7.5 


2660 


0.049 


16.473 


2955 


0.081 


16.456 


15 


Oil 


8 


38:06.120 


-53:38:10.55 


9.0 


M7.0 


2720 


0.050 


16.651 


2925 


0.076 


16.576 


15 


041 


8 


36:54.383 


-53:45:42.10 


35.2 


M5.5 


2925 


0.073 


14.988 


3266 


0.179 


14.953 


20 


001 


8 


40:34.407 


-52:30:38.65 


32.0 


M5.0 


3010 


0.089 


16.681 


2947 


0.079 


16.489 


20 


009 


8 


38:36.649 


-52:27:47.16 


20.2 


M5.0 


3010 


0.089 


16.873 


2911 


0.075 


16.631 


20 


012 


8 


37:38.539 


-52:29:35.37 


9.2 


M4.5 


3095 


0.110 


16.645 


2956 


0.081 


16.453 


20 


014 


8 


37:34.755 


-52:27:02.90 


27.6 


M5.0 


3010 


0.089 


16.307 


3026 


0.093 


16.161 


20 


018 


8 


37:58.411 


-52:20:30.69 


22.0 


M5.0 


3010 


0.089 


16.907 


2890 


0.072 


16.710 


20 


022 


8 


38:47.282 


-52:44:32.61 


17.0 


M6.0 


2840 


0.065 


16.662 


2927 


0.077 


16.567 


20 


023 


8 


39:22.724 


-52:50:34.42 


20.0 


M2.0 


3510 


0.413 


16.433 


2998 


0.087 


16.280 


20 


024 


8 


39:33.260 


-52:47:10.40 


21.8 


M6.5 


2780 


0.056 


16.520 


2958 


0.081 


16.445 


20 


028 


8 


40:15.153 


-52:40:24.56 


26.5 


M7.5 


2660 


0.049 


16.494 


2955 


0.080 


16.459 


20 


029 


8 


40:16.671 


-52:36:58.32 


21.4 


M6.5 


2780 


0.056 


17.571 


2740 


0.055 


17.329 


20 


033 


8 


39:27.836 


-52:32:58.62 


13.1 


M6.5 


2780 


0.056 


14.856 


3297 


0.197 


14.783 



Note. — The error on the determination of masses and effective temperature based on spectroscopy are the following : 
ATgff = 190 K and AM = 0.03 Mq for spectra with SNR > 10, AT^g = 320 K and AM = 0.04 Mq for spectra with 5 < SNR < 10. 
This would corresponds to an error on the spectral determination of 1 and 1.5 for spectra with SNR > 10 and 5 < SNR < 10 
respectively. 



Table 5. Object from iBarrado y Navascues et al.l (120041 ) which we also observed in our spectroscopic follow-up 



Field 


ID 


RA 


DEC 


SNR 


SpT 


Teff M 


NAME 


SpT Teff 


20 


022 


8:38:47.282 


-52:44:32.61 


17.0 


M6.0 


2840 0.065 


CTIO-062 


M6.0 2800 


Note. 




The SpT and 


Tefi given 


in the 


last two 


columns arc 


based on the 


spectra from 



ISarrado v Navascues et all ||2004|) 



Table 6. Spectroscopic data, photometric and spectroscopic membership status 



Field 


ID 


SNR 


SpT 


EWiKa) 


EW{Na.l 8182 A) 


EW{T<!aI 8194 A) 


RV (km/s) 


EW{Lil) 


Spec member ? 


1 c; 




D. i 


i\/ro n 
iViz.U 




1 01 








IN W 


1 c; 

io 


UU4 


on v 


iVio.O 


1 .OO 


Q 1 o 








1\TM 

IN yj 


io 


UUo 


Q Q 
O.O 


iVi / . 


4.oU 










Y rjo 


io 


m n 
UiU 


o.U 


IVi i . O 




U. ( o 








L\ yj 


io 


ni 1 
Uil 


Q n 
y.u 


IVl 1 .U 


-o.cSU 


O 1 
Z. iO 


i. 1 ( 








15 


041 


35.2 


M5.5 


12.66 


1.02 


1.75 






NO 


20 


001 


32.0 


M5.0 




1.21 


2.88 


15.11( 8.39 ) 


0.8 


YES 


20 


009 


20.2 


M5.0 


-17.56 


2.65 


2.19 


-18.28(12.37 ) 


<0.1 


YES 


20 


012 


9.2 


M4.5 


15.18 






11.31(27.70 ) 




YES 


20 


014 


27.6 


M5.0 


13.24 






23.33( 6.88 ) 


<0.1 


YES 


20 


018 


22.0 


M5.0 


21.36 


1.91 


3.32 


19.32( 7.78 ) 




YES 


20 


022 


17.0 


M6.0 


13.03 


1.55 


2.61 


14.95( 8.14 ) 




YES 


20 


023 


20.0 


M2.0 


-6.12 




1.72 


94.74( 6.90 ) 




NO 


20 


024 


21.8 


M6.5 




1.29 




27.91( 7.35 ) 




YES 


20 


028 


26.5 


M7.5 




1.23 




11.78( 7.48 ) 




NO 


20 


029 


21.4 


M6.5 




2.38 


2.33 


16.25( 6.59 ) 


1.3 


YES 


20 


033 


13.1 


M6.5 


5.80 




1.96 


6.23( 7.04 ) 




NO 



Note. — A value of <0.1 as equivalent width of Lithium indicates the presence of a feature at 6708 A but no reliable measurement 
of equivalent width can be done. 



Table 7. Stellar parameters obtained from spectra. 



Field 


ID 


RA 


DEC 


SNR 


Rc 


ilc-815/20 


SpT 


^eff 


M 


Teff (phot) 


M (phot) 


15 


005 


8:36:18.241 


-53:25:57.60 


8.3 


21.048 


3.245 


M7.5 


2660 


0.049 


2654 


0.048 


20 


029 


8:40:16.671 


-52:36:58.32 


21.4 


20.076 


2.505 


M6.5 


2780 


0.056 


2740 


0.055 



